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ABSTRACT 

We make black hole (BH) merger trees from Millennium and Millennium-II simula¬ 
tions to find under what conditions 10® Mq SMBH can form by redshift 2 : = 7. In order 
to exploit both: large box size in the Millennium simulation; and large mass resolution 
in the Millennium-II simulation, we develop a method to combine these two simu¬ 
lations together, and use the Millennium-II merger trees to predict the BH seeds to 
be used in the Millennium merger trees. We run multiple semi-analytical simulations 
where SMBHs grow through mergers and episodes of gas accretion triggered by major 
mergers. As a constraint, we use observed BH mass function at redshift z=6. We find 
that in the light of the recent observations of moderate super-Eddington accretion, 
low mass seeds (100 M©) could be the progenitors of high-redshift SMBHs {z ~ 7), 
as long as the accretion during the accretion episodes is moderately super-Eddington, 
where (/Edd = 3.7) is the effective Eddington ratio averaged over 50 Myr. 

Key words: stars: Population HI galaxies: high-redshift quasars: supermassive black 
holes dark matter 


1 INTRODUCTION 

Super-massive black holes (SMBHs) with masses 10® Mq 
to 10^® Mq populate centres of spiral galaxies with central 
bulge and centres of massive elliptical galaxies (Kormendy 
& Richstone 1995). These SMBHs grow through mergers 
and episodes of gas accretion. During accretion phase, black 
holes (BHs) are indirectly observed as active galactic nu¬ 
clei (AGNs) or quasars, where quasars powered by SMBHs 
with masses ~ 10® M© were observed as early as 2 ~ 6 
(e.g., Willott, McLure & Jarvis 2003; Fan 2006; De Rosa 
et al. 2014). Those observations showed that BHs had 
very little time (< IGyr) to grow from initial seeds (Fan 
et al. 2001). An extreme example is SMBH with mass 

2 X 10® M 0 detected at 2 = 7.085 (Mortlock et al. 2011). 
It would take a 10® Mq BH seed and continuous (many 
e-folding times) gas accretion at the Eddington limit, to 
grow such massive BH so early. Starting with lower mass 
BH seeds would require super-Eddington accretion. Obvi¬ 
ous question is: how massive were the BH seeds from which 
SMBHs formed so early and what are the key SMBH growth 
parameters? 

Several mechanisms for BH seeds formation have been 
proposed. We will briefly discuss the three most popular 
mechanisms. 

Pop HI stars, the first metal-free stars that started to 
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form at 2 ~ 20, might have left BH seeds that grew to 
SMBHs observed in the first billion years after the Big Bang 
(Madau & Rees 2001; Heger et al. 2003; Wise & Abel 2005). 
Masses of those BH seeds depend on the masses of Pop 
HI stars, which could be in range 60 — 300 Mq (Bromm 
et al. 2009), or even 1000 Mq (Hirano et al. 2014). Some 
simulations showed that Pop HI stars were not so massive 
and that they formed in binary systems or in clusters (Turk, 
Abel & O'Shea 2009; Clark et al. 2011; Stacy et al. 2012), 
but other authors showed that Pop HI stars were formed 
isolated (O'Shea & Norman 2008; Hirano et al. 2014). It is 
typically assumed in the literature that Pop HI stars form 
BH seeds with masses close to 100 Mq. 

Previous attempts to form ~ 10® Mq SMBHs at 2 ~ 6 
from 100 Mq seeds required continuous accretion close to 
or exceeding the Eddington limit and radiative efficiencies 
of e < 0.1 (Haiman & Loeb 2001; Tyler, Janus & Santos- 
Noble 2003; Volonteri & Begelman 2010; Whalen & Fryer 
2012; Johnson et al. 2012, 2013). Even with combination of 
BHs mergers and accretion, low-mass BH seeds must grow 
at the Eddington limit for a significant fraction of time or 
they need to have early stages of super-Eddington accre¬ 
tion in order to explain the observed high-redshift quasars 
(Yoo & Miralda-Escud 2004; Volonteri & Rees 2005, 2006; 
Li et al. 2007; Pelupessy, Di Matteo & Ciardi 2007; Sijacki, 
Springel & Haehnelt 2009; Tanaka & Haiman 2009; Tanaka, 
Perna & Haiman 2012; Madau, Haardt & Dotti 2014). There 
are theoretical uncertainties whether BH seeds can sustain 
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such high accretion rates for a long time (Milosavljevic, 
Couch & Bromm 2009, Alvarez, Wise & Abel 2009, Jeon 
et al. 2011). It is even more difficult to form ~ 10® Mq 
SMBH s at 2 ~ 7 considering that the time available for 
accretion is even smaller than at 2 ~ 6. 

Another possible mechanism for BH seed formation is 
direct collapse of gas into BH. This process produces BH 
seeds with masses in range ~ 10"^ — 10® M© (Loeb & Rasio 
1994; Eisenstein & Loeb 1995; Oh & Haiman 2002; Bromm 
& Loeb 2003; Koushiappas, Bullock & Dekel 2004; Begel- 
man, Volonteri & Rees 2006, Lodato & Natarajan 2006, 
Begelman, Rossi & Armitage 2008). It has been shown that 
massive BH seeds formed by direct collapse can reach IO^Mq 
at 2 > 6 (e.g., Dijkstra et al. 2008; Agarwal et al. 2012; 
Bonoli, Mayer & Callegari 2014; Petri, Ferrara & Salvaterra 
2012; Johnson et al. 2012, 2013). However, BH seeds in 
these models may still need to accrete at nearly the Ed¬ 
dington limit for a significant fraction of time. The main 
challenge in direct collapse models is to avoid fragmenta¬ 
tion and star formation. Fragmentation will not occur if gas 
has no mechanism for efficient cooling, that is if there are 
no metals and if formation of H 2 molecules is prevented via 
strong UV background. Only extremely rare haloes satisfy 
this condition. Direct collapse is expected to occur in haloes 
with virial temperature Tvir ^ 10‘^K and masses 10^ —IO^Mq 
at 10 < 2 < 20. Recently, models have been proposed where 
direct collapse might be possible even without a strong UV 
radiation background (Inayoshi & Omukai 2012; Tanaka & 
Li 2013), but Visbal, Haiman & Greg (2014) have shown 
that those models are not viable because molecular cooling 
will still occur, as the gas density increases, which leads to 
fragmentation. Alternatively, major mergers may lead to a 
rapid gas inflow. In such case, turbulence may be the in¬ 
hibitor of fragmentation, and the requirement of metal-free 
gas may be relaxed (Mayer et al. 2010). 

The third possible mechanism for BH seed formation 
is collapsing star clusters. Devecchi & Volonteri (2009) sug¬ 
gested model where mergers of Pop H stars lead to the for¬ 
mation of a very massive star which collapses into a BH with 
mass ~ 10® Mq, at 2 ~ 10 — 20. Davies, Miller & Bellovary 
(2011) proposed a model of a star cluster which contains only 
BHs and main sequence low mass stars. More massive BHs 
sink to the centre and form BH cluster which further col¬ 
lapses and forms massive BH with mass 10® Mq, at 2 > 10. 
These seeds can grow via mergers and accretion to form 
SMBHs. Some authors showed that such seeds could ex¬ 
plain observed quasars at 2 < 5 but have trouble explaining 
quasars at higher redshifts (Ebisuzaki et al. 2001; Volonteri, 
Haardt & Madau 2003; Islam, Taylor & Silk 2003). 

BH mergers could contribute to the SMBH growth. 
Johnson et al. (2013) argued that growth of the most mas¬ 
sive BH at high redshift is solely due to the accretion and 
that mergers can be neglected. Their conclusion is based on: 
the fact that only a small number of high-redshift haloes 
are capable of hosting Pop HI stars; Pop III star forma¬ 
tion rate is reduced due to Lyman-Werner (LW) radiation 
field; and the results of the recent large-scale cosmologi¬ 
cal simulations tracking the buildup of SMBH, where it is 
claimed that mergers do not affect the growth of SMBHs 
(DeGraf et al. 2012). Simulation by DeGraf et al. (2012) 
does show that negligible fraction of high-redshift SMBH 
mass is acquired in mergers but this is a consequence of 


the poor mass resolution. Their simulation does not form 
haloes smaller than 5 x 10®®/i~® M©. While it is true that 
Pop HI star formation can not take place in smallest haloes, 
it does occur in haloes with Tvir ~ 1000 — 2000K which 
corresponds to Mhaio ^ 10® M© at redshift 2 > 6. Also, 
haloes with Mhaio (J 10® M© are most likely self-shielded 
from LW radiation due to the large H 2 column density 
(Wise & Abel 2008). Hence, it is clear that haloes with mass 
10® M© < Mhaio ^ 10®® M© host Pop HI stars and if these 
stars produce massive BHs, then their mergers should be 
relevant to the growth of high redshift SMBHs. 

We investigate if light BH seeds (100 M©) planted into 
haloes of Millennium simulation (Mhaio > 10®® M©) and 

Millennium-II simulation (Mhaio > 10® M©) can grow into 
SMBHs that have been observed at 2 ~ 7. 

We also investigate if gravitational wave recoil could 
prevent the formation of SMBH. During a BH merger grav¬ 
itational wave radiation is produced. Asymmetric emission 
of gravitational radiation can lead to BH kick. Gravitational 
waves carry a non-zero net linear momentum, which es¬ 
tablishes a preferential direction for the propagation of the 
waves and the centre of mass of the binary recoils in the op¬ 
posite direction (Redmount & Rees 1989). The magnitude 
of the gravitational wave recoil depends on the mass ratio 
of BHs, the spin magnitude and orientation with respect to 
the binary orbital plane and the eccentricity of the orbit 
(Campanelli, Lousto & Zlochower 2007; Schnittman & Buo- 
nanno 2007; Baker et al. 2002). Gravitational wave recoil 
can eject a newly formed BH from the host halo if the BH 
speed is larger than the escape velocity from the halo centre. 
BHs can be kicked with a speed as large as ~ 4000 km s“® in 
special orbital configurations (Herrmann et al. 2007; Gonza¬ 
les et al. 2007a, 2007b; Campanelli et al. 2007; Schnittman 
& Buonanno 2007; Koppitz et al. 2007). At high redshift, 
dark matter haloes generally have small masses and thus 
small escape velocities, so BH with speed 150 km s“® 
can be ejected even from the most massive haloes at red¬ 
shift 2 is 11 (Merritt et al. 2004; Micic, Abel & Sigurdsson 
2006; Volonteri 2007; Schnittman 2007; Sesana 2007; Volon¬ 
teri, Gultekin & Dotti 2010; Micic, Holley-Bockelmann & 
Sigurdsson 2011). This effect may play a major role in sup¬ 
pressing the growth of SMBH through mergers. 

1.1 Growth Parameters and Their Values 

Large scale structure formation and galaxy dynamics lead 
to BH mergers. If those mergers are ignored, BH growth de¬ 
pends on three gas accretion parameters. Those parameters 
are radiative efficiency, Eddington ratio and the time that a 
BH spends accreting. 

Radiative luminosity of a BH, L, which is accreting at 
a rate Mbh, is given by T = e Mbh c ®/(1 — e), where c is the 
speed of light and e is the radiative efficiency. Eddington ra¬ 
tio is /Edd = where LEdd = 1-26 x 10®®(^^^)[erg s“®] 

•'-'Edd 

is Eddington luminosity. The accretion rate at which a BH 
will radiate at a given Eddington ratio is given by: 

Mbh = ( 1 ) 

After integration the final BH mass Mbh, as a function of 
its initial mass Mbh.o is: 
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Mbh = Mbh.o X exp [ ^ ( 2 ) 

e iBdd 

where ^Bdd = 450 Myr, if i if are the ages of the universe 
when the BH attains its final mass and at the time of seed 
formation, respectively (Johnson et al. 2013). 

Radiative efficiency is the efficiency of conversion of 
rest-mass into energy during accretion and it depends on 
the BH spin. Radiative efficiency can take values from 0.057 
for accretion on to Schwarzschild BHs to 0.42 for fast ro¬ 
tating Kerr BHs (Shapiro 2005). Mean value of radiative 
efficiency for quasars can be estimated by comparing the lo¬ 
cal SMBH mass density with the total AGN luminosity per 
unit volume in the Universe integrated over time (Soltan 
1982). This is Soltan’s argument. Previous works based on 
Soltan’s argument showed that the mean value of radia¬ 
tive efficiency is e is 0.1 (Elvis, Risaliti & Zamorani 2002; 
Yu & Tremaine 2002; Davis & Laor 2011). Some authors 
have found that radiative efficiency changes with the red- 
shift (Wang et al. 2009; Li, Wang & Ho 2012) and increases 
with the mass of the accreting BH (Davis & Laor 2011; 
Shankar, Weinberg & Miralda-Escud 2011; Li, Wang & Ho 
2012 ). 

It is usually assumed that BH luminosity during ac¬ 
cretion cannot be greater than the Eddington luminosity. 
However, evidence for super-Eddington accretion has been 
growing recently. Kelly & Shen (2013) used a sample of ~ 
58,000 quasars at s ~ 0.3 — 5 from Sloan Digital Sky Sur¬ 
vey (SDSS) DR7 catalogue (Schneider et al. 2010) to esti¬ 
mate their Eddington ratios. They found that the highest 
observed Eddington ratio for quasars is /Bdd ~ 3, but that 
those quasars which radiate above the Eddington limit are 
rare. Du et al. (2014) observed three quasars, Mrk 335, 
Mrk 142 and IRAS F12397+3333, using the 2.4-m Shangri- 
La telescope at the Yunnan Observatory in China. One 
of their goals was to measure BH masses and Eddington 
ratios. They found that the lower limits on the Edding¬ 
ton ratios for these objects are 0.6, 2.3, and 4.6. Page et 
al. (2014) have shown that soft X-ray spectrum of the high¬ 
est redshift quasar yet found, ULAS J112001.48-1-064124.3 
at « = 7.085 (Mortlock et al. 2011), obtained with Chandra 
and XMM — Newton, suggests that the quasar is accreting 
above the Eddington limit, /Bdd = 5 I 4 ®. Their hndings of 
moderate super-Eddington accretion are consistent with the 
Eddington ratio estimated from the UV luminosity of that 
object, /Bdd = l.Slg ® (Mortlock et al. 2011). In order to 
fit the observed statistics of far-infrared and X-ray spectra 
of AGNs at 2 > 2 Lapi et al. (2006, 2014) assumed model 
where Eddington ratio depends on the redshift. During the 
exponential growth of the BH, the maximum Eddington ra¬ 
tio is /Bdd ~ 4 for 2 = 6 , and /Bdd ~ 1 for 2 = 2, with 
constant radiative efficiency e = 0.15. Similar results have 
been found by Li (2012), who has also explored the possibil¬ 
ity of having a short super-Eddington accretion followed by 
a sub-Eddington accretion in order to explain the presence 
of BHs with masses ~ 10® M© SMBHs at 2 ~ 6 . 

A recent analysis of BH scaling relations showed that 
the normalization of the BH mass-bulge relation should be 
increased by a factor of 5, from the previously accepted 
value of Mbh = 0.1%Mbuige to Mbh = 0.5%Mbuige. This 
increases the local mass density in BHs by the same fac¬ 


tor and decreases the required mean radiative efficiency to 
values that cannot be reasonably explained in terms of lu¬ 
minous thin-disc accretion. This may be evidence for the ra- 
diatively inefficient super-Eddington accretion (Novak 2013; 
Kormendy & Ho 2013). These works have showed that 
quasars, at least at some point in their evolution, can ac¬ 
crete at the super-Eddington luminosities. Recent numerical 
simulations (McKinney et al. 2014; Sadowski et al. 2013) 
of super-Eddington accretion suggest that those sources are 
most likely to be characterized by strongly collimated out¬ 
flows or jets. Accretion in these simulations is mildly super- 
Eddington, /Edd = 1—10. 

Theoretical works by Volonteri & Rees (2005, 2006) as¬ 
sumed early stages of super-Eddington quasi-spherical ac¬ 
cretion estimated using the Bondi-Hoyle formula (Bondi & 
Hoyle 1944). In this case, when the inflow rate is super¬ 
critical, the radiative efficiency drops. Hence, the Eddington 
luminosity is not greatly exceeded (I/Bdd = eMnddC®). The 
accretion rate is initially super-critical by a factor of 10 and 
then grows up to a factor of about 10"^. They showed that 
BH seeds with initial mass of 1000 Mq could explain SMBH 
at 2 ~ 6 if they go through early phases of super-Eddington 
accretion and BH mergers. Similar approach has been re¬ 
cently argued by Volonteri & Silk (2015). They showed that 
short-lived intermittent episodes of super-Eddington accre¬ 
tion (/Edd > 10) may increase the BH mass by several orders 
of magnitude in ~ 10^ years. The authors assumed slim disc 
solution where the luminosity during accretion depends log¬ 
arithmically on the accretion rate. Since (/nud = ^ ); 

if the effective radiative efficiency is low, accretion rate 
can be highly super-Eddington while the emergent lumi¬ 
nosity is only mildly super-Eddington and feedback is lim¬ 
ited. Madau, Haardt & Dotti (2014) extended these works 
and showed that light BH seeds (100 Mq) could explain 
quasars with 10® Mq SMBH observed at 2 ^ 6 if they have a 
few episodes of super-Eddington accretion via slim accretion 
disc. Eddington ratio in these works is ~ 1 — 10. 

Dehnen & King (2013) proposed different mechanism 
for super-Eddington accretion. Momentum driven feedback 
from an accreting BH gives significant orbital energy but 
little angular momentum to the surrounding gas. Once cen¬ 
tral accretion drops, the feedback weakens and gas falls back 
towards the BH, forming a small scale accretion disc. The 
feeding rates into the disc typically exceed Eddington by 
factor of a few. 

Based on the mentioned works it can be concluded that 
super-Eddington accretion is inevitable for the growth of 
SMBHs. 

Another important question is for how long accretion 
can be sustained? Quasars lifetimes can be estimated us¬ 
ing Soltan’s argument. Amount of matter accreted on to 
quasars during their lifetime, represented by the luminos¬ 
ity density due to accretion, should be less than or equal 
to the space density of remnant BHs in the local Universe 
(Soltan 1982). This method depends on the value of radia¬ 
tive efficiency. Another way to estimate this parameter is to 
ask what value of quasars lifetime is required if all bright 
galaxies go through a quasar phase. This approach is not 
sensitive to the value of radiative efficiency, but it is af¬ 
fected by the assumed value of Eddington ratio and galaxy 
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mergers (Martini 2004). Computed values for the duration 
of accretion using this approach are from 10 ® to 10 ® yr. 

Typical accretion time for majority of quasars is 
Salpeter’s time (Salpeter 1964). Salpeter’s time is e-folding 
time-scale for SMBH growth and its typical value is 
~ 50 Myr: 

ts = M/M = 4.5x 10^(^)(^)-\ (3) 

A single massive BH seed (10® — 10® M©) needs to con¬ 
tinuously accrete at Eddington luminosity for 380-500 Myr 
(which corresponds to 9 - 11 e-folding times) to be able to 
grow to > 2 X 10® Mq at redshift 2 6. Light BH seed 

(100 Mq) requires 19 e-folding times or 840 Myr of accre¬ 
tion. Alternatively, it would take 280 Myr at Eddington ra¬ 
tio 3 (also 19 e-folding times). This means that it would be 
very hard to grow SMBH from light seeds at the Eddington 
limit by 2 = 6 and impossible by 2 = 7. Some authors as¬ 
sumed such long quasar lifetimes in order to match observed 
quasars number density at high redshifts (Haiman & Loeb 
2001; Tyler et al. 2003; Sijacki et al. 2009; Tanaka & Haiman 
2009; Tanaka et al. 2012; Johnson et al. 2012, 2013). 

Continuous accretion at or above the Eddington limit 
might be problematic due to the feedback. Radiation and 
kinetic power in matter outflows could be so strong that 
once the accreting BH gets too big, it blows out all of the 
gas in the centre of its host galaxy, shutting down the accre¬ 
tion on to the BH (Coppi 2003). If continuous growth at the 
Eddington limit cannot be maintained, the main alternative 
to these models are short periods of super-Eddington accre¬ 
tion (Volonteri & Rees 2005, 2006; Dehnen & King 2013; 
Volonteri & Silk 2015; Madau, Haardt & Dotti 2014). How¬ 
ever, observations show that, in most cases, quasars lifetimes 
are comparable with Salpeter’s time and Eddington ratios 
are between 0.1 and 1. Time that SMBH spends accret¬ 
ing depends on radiative efficiency, Eddington ratio and BH 
mass (Martini 2004). Yu & Tremaine (2002) showed that 
quasar lifetime is a function of BH mass. They found that 
the mean lifetime is 3 — 13 x 10^ yr for e = 0.1 — 0.3 and 
10® < Mbh < 10® Mq. It is obvious that if one assumes 
these typical values of the growth parameters, it would be 
impossible to grow SMBH at high redshift even if we start 
with seeds as massive as 10^ M©. Accretion would have 
to be either super-Eddington or prolonged for more than 
Salpeter’s time (many e-folding times instead of one). 

Our goal is to examine if BH mergers combined with 
accretion episodes in merger trees of Millennium simulation 
(Springel et al. 2005) and Millennium-II simulation (Boylan- 
Kolchin et al. 2009) contribute to the growth of 10® M© 
SMBH at 2 ~ 7. If they do, then values of Eddington ratios 
and numbers of e-folding times during the merger induced 
accretion episodes in merger tree branches can be lowered 
to a more reasonable values. Growth of a SMBH in halo 
that undergoes several major mergers does not require con¬ 
tinuous accretion for many e-folding times. Since the rate 
of major mergers increases with redshift, SMBHs in high- 
redshift quasars can grow nearly continuously in a sequence 
of major mergers (Li et al. 2012, Tanaka 2014). Every ma¬ 
jor merger will reset previous accretion episode and trigger a 
new one which effectively reduces number of e-folding times. 


That in turn could increase the likelihood of low-mass BHs 
as SMBH seeds. 

In Section 2 we describe the method. In Section 3 we 
present our results. We summarize and discuss our results 
in Section 4. 


2 METHOD 

Our goal is to produce SMBH with mass > 10® M© at 
2 ~ 7 from light BH seeds using publicly available data 
from both Millennium simulation (Springel et al. 2005) and 
Millennium-II simulation (Boylan-Kolchin et al. 2009). As 
a constraint we use comparison between the observed and 
calculated BH mass function from our model at 2 ~ 6 . 


2.1 Millennium simulation and Millennium-II 
simulation 

Millennium simulation (Springel et al. 2005) is a large N- 
body simulation which follows 2163® particles within a peri¬ 
odic simulation cube of side length L = 500/i“® Mpc. Each 
simulation particle has mass 8.61 x 10® M©. The ACDM 
cosmology used for the Millennium simulation is: 

fitot = 1.0, Ilm = 0.25, fib = 0.045, flA = 0.75, h = 0.73, 
as = 0.9, Us = 1, 

where h is the Hubble constant at 2 = 0 in units of 
100 km s“® Mpc“^, as is the rms amplitude of linear mass 
fluctuations in 8h~^ Mpc spheres at 2 = 0, and Us is the 
spectral index of the primordial power spectrum. 

Millennium-II simulation (Boylan-Kolchin et al. 2009) 
uses the same cosmology and the same number of particles as 
the Millennium simulation, but in a five times smaller box 
(L = 100h“^ Mpc) and thus with 125 times better mass 
resolution. Each simulation particle in Millennium-II simu¬ 
lation has mass 6.885 x 10® M©. Millennium-II simulation 
uses GADGET-3 code, updated version of GADGET code 
(Springel, Yoshida & White 2001b, Springel 2005). 

Current dominant cosmological paradigm is CDM, 
which predicts bottom-up mode of structure formation, that 
is, small dark matter haloes forming first and then merging 
into larger haloes later in the fife of the Universe. 

Since Millennium-II simulation has 125 times better 
mass resolution than Millennium simulation, it allows as to 
track BH growth by mergers of highest redshift low mass 
haloes which are too small to be resolved due to the lower 
mass resolution of Millennium simulation. On the other 
hand, it is questionable if it is possible to produce a 10® M© 
SMBH at 2 ~ 7 in a simulation box with sides of length 
that Millennium-II simulation has due to low population of 
highest mass haloes. In order to have both small haloes and 
large box, we develop a method to use BHs produced in 
merger trees of Millennium-II simulation, as BH seeds in 
merger trees of Millennium simulation, which enable us to 
cover both haloes that are forming very early and abundant 
highest mass haloes later. 
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2.2 Combining Millennium-II and Millennium 
merger trees 

With 100 M© BH seeds placed in Millennium-II haloes we 
produce BH growth history. From it, because of the higher 
mass resolution of Millennium-II simulation, we have infor¬ 
mation about earliest halo formation history. 

First, we make merger trees which track dark matter 
halo merger history in Millennium-II simulation from red- 
shift a = 23.79 to 2 = 6.2. Halo is defined as self-bound 
structures with at least 20 particles. Merger history of BHs 
corresponds to the merger history of haloes assuming that 
each halo hosts one BH and that two BHs merge right after 
their host halo merge. We have distinguished between mi¬ 
nor and major mergers. Merger is major if ^ 0.3 for 

-^^halo,2 

Mhaio.i < A7haio.2. lu the case of a minor merger, mass of 
the newly formed BH is a simple sum of the previous BH 
masses. In the case of a major merger, newly formed BH is 
also accreting gas according to the accretion recipe described 
in following subsections. In each halo of Millennium-II sim¬ 
ulation we place one BH seed with 100 M© and assume fixed 
values for radiative efficiency and accretion time-scale in ev¬ 
ery accretion episode. Hence, the only variable in our model 
is the Eddington ratio. 

Next, we use Millennium simulation which has simula¬ 
tion box five times larger {L = 500/i“^ Mpc) in order to 
compare final BH masses in merger trees of this simulation 
to the observed BH mass function. We apply the same BH 
growth recipe to the haloes of Millennium simulation. The 
difference is, time the seed BH mass is the most common BH 
mass estimated in the merger trees of Millennium-II simu¬ 
lation. 

We use Millennium-II simulation to determine a typ¬ 
ical BH mass for a specific host halo mass. We do it by 
binning masses of resulting BHs at each snapshot and halo 
mass interval and then choosing central value of the most 
occupied bin. Hence, at every redshift we know the range 
of BH masses hosted by haloes with a specific mass. We se¬ 
lect the typical (most common mass) BH and use it as a 
BH seed in each newly formed halo of Millennium simula¬ 
tion at the matching redshift and halo mass obtained from 
Millennium-II, thus extending our sample to include highest 
mass haloes. Note that this is a conservative way to popu¬ 
late haloes since we are losing instances when BH mass is in 
reality more massive than the common one. 

2.3 Merger tree 

First, we use Millennium-II data base to select haloes with 
masses > 10^° M© at redshift 2 = 6.2. For all selected haloes 
we find their progenitor haloes, i.e. haloes that have merged 
at the previous snapshot to form selected haloes. Haloes have 
merged if they have the same descendant halo. We repeat 
this procedure up to redshift 2 = 23.79 where first mergers 
are recorded. Once we find all of the haloes that take part 
in formation of the selected haloes, we form merger trees for 
each of them. 

Final haloes at redshift 2 = 6.2 are called main haloes. 
Mergers of the main halo form the main branch of the merger 
tree. Other haloes that have merged with the main halo are 
side haloes and their previous mergers form side branches. 
Every merger tree contains main branch and side branches. 


We follow not only mergers of the main halo, but also take 
into account mergers in all side branches. 

In a single merger event the most massive progenitor 
halo is called primary halo and other progenitor haloes are 
satellite haloes. In the Millennium-II simulation output is 
printed with time resolution Az ~ 1, so the common case 
is that large number of haloes have the same descendant 
halo (case with largest number of progenitors being 186). 
The question is which halo is the primary halo and which 
haloes are satellite haloes? To solve this we first find the 
most massive progenitor halo (primary halo) and then using 
the comoving coordinates of the centre of mass, we calculate 
distances to other progenitor haloes (satellite haloes), sort 
them in the order from closest to furthest away and assume 
that they have merged in that order. We calculate halo mass 
ratio of the merging haloes and distinguish minor and major 
mergers. 

We repeat the same procedure for haloes in Millennium 
simulation that have masses > 10^^ M© at redshift z = 6.2. 

2.4 Parameters choice and initial BH masses 

Final BH mass depends on the initial BH mass, Eddington 
ratio, radiative efficiency and accretion time-scale (equation 
2.). In our model, every major merger leads first to the for¬ 
mation of a new BH after which gas accretion phase is trig¬ 
gered. 

Since we are only interested in the initial and the fi¬ 
nal BH mass before and after every accretion episode, we 
do not strictly insist on any particular accretion model. The 
accretion parameters in our model should be regarded as av¬ 
eraged over the accretion time-scale no matter what the ac¬ 
tual BH growth model is. These averaged parameters could 
be the real accretion parameters in the classical ’thin disc’ 
accretion model, or they could be an average of a sequence 
of short-lived intermittent phases of super-Eddington accre¬ 
tion (with Eddington ratios greater than the averaged one) 
in the ’slim disc’ accretion model (Volonteri & Silk 2015). 
That is why we refer to Eddington ratio and radiative ef¬ 
ficiency as effective Eddington ratio and effective radiative 
efficiency in our model. 

For effective radiative efficiency we choose commonly 
accepted value e = 0.1 (Elvis, Risaliti & Zamorani 2002; 
Yu & Tremaine 2002; Davis & Laor 2011). Every accretion 
episode is limited to 50 Myr, which is ~ Salpeter’s time for 
accretion at the Eddington ratio f^dd = 1- The only free pa¬ 
rameter in our model is the effective Eddington ratio. We as¬ 
sign a fixed value of this parameter to each accretion episode 
in one simulation run. Accretion episode can be shorter than 
50 Myr in two cases. First, if a new major merger occurs 
before that time has passed, accretion is reset and new ac¬ 
cretion episode begins. Secondly, if BH mass exceeds 0.08 
per cent of the host halo mass, accretion will be stopped 
and cannot be triggered again. This constraint comes from 
BH-dark matter halo mass scaling relation in the local Uni¬ 
verse. Mass in baryons is approximately 16 percent of the 
dark matter mass. If a typical BH mass to bulge mass ratio 
is 0.5 per cent (Kormendy & Ho 2013, Novak 2013), then 
8 X 10“^ of the halo mass is the gas that BH can accrete. 

In each newly formed halo of Millennium-II simulation 
we place one BH with initial mass of 100 M© and run one 
simulation with fixed accretion parameters. We use haloes of 
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Millennium-II simulation to estimate what are the common 
masses of BHs that populate haloes of a certain mass at each 
redshift. Then we use those common BHs as seeds in newly 
formed haloes of Millennium simulation and apply the same 
BH growth recipe with the same accretion parameters as we 
have used in Millennium-II simulation. 

We use the observed BH mass function at z ~ 6 as a 
constraint for our model. We run a set of semi-analytical 
simulations for different values of the effective Eddington 
ratio. In each run we assign the same values for the effective 
Eddington ratio to each accretion episode. We make sure 
not to overproduce SMBH at 2 = 6.2 (Millennium snapshot 
closest to redshift 2 = 6 ) and once this condition is satisfied, 
we check if we have 10® Mq SMBH at 2 = 7 for the specific 
choice of the effective Eddington ratio. 


3 RESULTS 

We find that if the effective Eddington ratio is /sdd = 3.7 
both conditions of our model are satisfied: BH mass function 
is consistent with the observed BH mass function at 2 ~ 6 
(Willott et al. 2010b) and our merger tree produces 10® Mq 
SMBH at 2 = 7. 

Fig. 1 shows mass function of BHs that populate > 
10^^ Mq haloes of Millennium Simulation at 2 = 6.2. BH 
masses are the result of one semi-analytic simulation in 
which we chose initial BH masses to be 100 Mq, effec¬ 
tive radiative efficiency e = 0.1, effective Eddington ratio 
/bjj = 3.7 and each accretion episode is limited to 50 Myr. 
In order to calculate BH mass function in our model we 
bin masses of all BHs in haloes of Millennium simulation 
at redshift z=6.2. Points depicted by plus symbols in Fig. 1 
represent numbers of BHs in each bin per Mpc®. Bins have 
width of 0.1 dex (in logarithmic scale). We compare our BH 
mass function to the BH mass function at 2 = 6 given by 
Willott et al. (2010b) (dashed and solid lines). Willott et al. 
(2010b) modelled BH mass functions to produce luminosity 
functions which are then fitted to the observed luminosity 
function of Canada-France High-z Quasar Survey (CFHQS) 
and SDSS quasars at 5.74 < 2 < 6.42 (Willott et al. 2010a). 
Dashed lines in Fig. 1 represent their upper and lower limits, 
while the solid line is their best fitting to the data. 

Kelly & Shen (2013) showed that accretion with similar 
Eddington ratios exists in high redshift quasars. They found 
that the maximum observed value of Eddington ratio for the 
observed quasars is /Edd ~ 3. However, other authors sug¬ 
gest that even larger values of Eddington ratio are possible 
(/Edd = 4.6 (Du et al. 2014), or even /sdd = 10 (Collin & 
Kawagachi 2004)). 

Assuming that a moderate super-Eddington accretion 
is possible for a prolonged period of time (more than one 
e-folding time), in the further analysis we focus on the main 
halo with the most massive SMBH at redshift 2 = 7 and its 
progenitors. 

Fig. 2 shows the merger tree of a 10® Mq SMBH at 
redshift 2 = 7 and effective /Edd = 3.7. It follows the en¬ 
tire growth of a SMBH through all mergers and accretion 
episodes, as a function of the age of the Universe. Black cir¬ 
cles represent BH masses in side haloes (side branches) at 
the snapshot-times, while red asterisks represent BH masses 
in the main halo (main branch). Dotted blue lines follow BH 



Figure 1. Mass function of BHs in our model (blue plus sym¬ 
bols), compared to BH mass function given by Willott et al. 
(2010b). Dashed lines show upper and lower limit from Willott et 
al. (2010b) while the solid line shows their best-fitting. 



_^^^^^^^^_I_^^^^_ 
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Figure 2. Merger tree of a 10® Mq SMBH at redshift 2 = 7 for 
/Edd = 3.7. It follows the entire growth of a SMBH through all 
mergers and accretion episodes, as a function of the age of the 
Universe. Black circles represent BH masses in side haloes at the 
snapshot-times, while red asterisks represent BH masses in the 
main halo. Dotted blue lines follow BH growth by minor mergers 
and solid red lines show major mergers which are always followed 
by the gas accretion. 

growth by minor mergers and solid red lines show growth 
through the gas accretion which occurs after every major 
merger. At the end, there is only one SMBH left. BH gains 
most of its mass in major mergers when accretion is trig¬ 
gered. 

Accretion episodes at t < 300Myr occur in haloes 
of Millennium-II simulation, while accretion episodes at 
t > 600Myr represent major mergers of haloes in Millen¬ 
nium simulation. The absence of accretion episodes between 
t = 300Myr and t = 600Myr is the consequence of the 
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Figure 3. Growth of the SMBH in the main halo by minor merg¬ 
ers (black solid line), major mergers (blue dashed line) and direct 
accretion (red dash-dotted line). Below 10^ Mq major mergers 
and accretion are equally important while above that mass SMBH 
growth is due to direct accretion. 

fact that we do not know the exact merger history of low 
mass haloes of Millennium simulation. We approximate BH 
masses in those haloes with typical BHs that populate haloes 
of Millennium-II simulation at the matching redshift. 

Each accretion episode lasts 50 Myr. Salpeter’s time for 
a BH accreting at the effective Eddington ratio /Bdd = 3.7 
is ~ 12Myr (equation 3), which means that e-folding time 
in each accretion episode is ~ 4. 

Super-Eddington accretion for such a long period of 
time might seem unrealistic, but it could be more possi¬ 
ble than the accretion at the Eddington limit for almost a 
billion years. To explain the formation of > 10® Mq SMBH 
with a former classical values of the accretion parameters 
(Mbh.o = 10 ® M 0 , e = 0.1 and f^dd = 1 ) it can be calcu¬ 
lated that continuous accretion for 500 Myr is needed (equa¬ 
tion 2.) and e-folding time is ~ 11. If we start with those 
parameters and smallest BH seed of 100 Mq, we calculate 
that e-folding time is ~19 and the accretion lasts for 840 
Myr in case where /sdd = 1 and 280 Myr in case where 
/sdd = 3. 

In our model, growing SMBH a.t z = 7 does not require 
e-folding time larger than 4 for 100 Mq BH seeds. This is a 
consequence of accretion being restarted in major mergers. 
Instead of one BH constantly accreting we have several BHs 
in shorter accretion episodes, where each of them is triggered 
in a major merger. SMBH in the main halo grows through 
two accretion episodes while the rest of the accretion occurs 
in side haloes (Fig. 2). Above 10^ Mq two accretion episodes 
occur in a side halo. BH in that side halo grows parallel 
with BH in the main halo. When the side and the main 
haloes merge side halo hosts more massive BH than the main 
halo because it has greater number of major mergers in its 
history. This approach increases the impact of mergers, and 
reduces importance of accretion, which in turn alleviates 
need for super-Eddington accretion with the large number 
of e-folding times. 

Fig. 3 represents accumulative contribution to the 
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Figure 4. Main and side haloes masses as a function of BHs 
mass ratio in both major and minor halo mergers. Red diamonds 
represent haloes where the sum of two merging BHs is > 10^ Mq. 
Those mergers are the most important in BH growth and they 
are not affected by gravitational wave radiation. Mergers of BHs 
with mass ratio ~ 1 might suppress the growth of SMBH by 
gravitational wave recoil since they reside in low-mass haloes. 
BHs that form in larger haloes are protected from ejection by the 
large gravitational potentials. 


SMBH growth from minor mergers (black solid line), major 
mergers (blue dashed line) and direct accretion (red dash- 
dotted line). It follows the growth of SMBH in the main halo. 
Major mergers have significant contribution below 10^ Mq 
because BHs in side haloes can grow to high masses before 
they merge with the BH in the main halo. Above 10^ Mq 
SMBH growth is dominated by direct accretion. 

We also investigate if gravitational wave recoil could 
prevent the formation of SMBH. Fig. 4 shows the mass ratio 
of the merging BHs in the host haloes oi z = 7 merger tree. 
It shows the main and side haloes masses as a function of 
BH mass ratio. Our model might be sensitive to the choice of 
kick velocity. Most of the BH mergers with mass ratios close 
to one occur in low mass haloes. As low mass haloes have 
small gravitational potentials, hnal BH after the merger can 
easily be ejected. If BH survives these mergers and settles 
in a larger halo, then the larger gravitational potential will 
protect it from any following merger. Red diamonds repre¬ 
sent haloes where the sum of two merging BHs is > 10^ Mq. 
Those mergers have the largest influence on the BH growth 
and since they occur in large haloes and have low mass ratio, 
they are not affected by the gravitational wave recoil. Fig. 
4 also shows that mass ratio of merging BHs decreases in 
larger mass haloes. 

To further investigate the influence of the gravitational 
wave recoil, we calculate and compare kick value with the 
value of the halo escape velocity. 

Kick velocities are taken from Micic et al. (2011). The 
authors used parametrized fit of Campanelli, Lousto & Zlo- 
chower (2007) to calculate the kick velocity as a function of 
the merging BHs mass ratio, BHs spin and the alignment 
to the orbital angular momentum. Kick velocity is (their 
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Figure 5. Number of BH mergers as a function of BH mass. Grey 
filled histogram shows BHs whose kick velocities are larger than 
the escape velocities from their host halo, so they are sensitive 
to the gravitational wave recoil. Blue histogram with horizontal 
lines shows BHs that are not affected by gravitational wave recoil. 


equations 7-10): 


Wkick = [(Wm + v±cos^)^ + {v± sin^)^ + 


where 




(1 - g) 
(i + g)® 


l + B 


g 

(i + gf. 


(4) 

(5) 


with NFW profile (Navarro, Frenk & White 1996) with a 
concentration parameter c = 4 out to the virial radius, as 
expected for a newly formed dark-matter halo (Wechsler et 
al. 2002). Under those assumptions the escape velocity from 
the halo’s centre is Vesc ~ 2.8vc. 

Fig. 5 shows the number of BH mergers as a function 
of the mass of the newly formed BH. If BH kick veloc¬ 
ity is larger than the escape velocity from the host halo, 
newly formed BH could be ejected from the halo (grey filled 
histogram). Blue histograms with horizontal lines represent 
BHs whose kick velocities are smaller than the escape veloc¬ 
ities. The figure shows that low-mass BHs reside in low mass 
haloes and have their kick velocities larger than the escape 
velocities from their host haloes, so they are sensitive to the 
gravitational wave recoil. On the other hand, we use simple 
model in which we assume that each newly formed halo in 
Millennium-H simulation hosts a BH with the initial mass 
of 100 M 0 . Possibility that two merging BHs have exactly 
the same mass is very small. First BHs have initial mass 
function (IMF), with masses in wide range. IMF for Pop HI 
BH remnants is not fully understood due to uncertainty in 
primordial gas fragmentation during Pop HI star formation. 
Pop HI stars with masses in range ~ 25 — 140 Mq are be¬ 
lieved to leave BH remnants with masses Mbh ~ 10 — SOMq 
while more massive stars, > 260 Mq, leave BHs with masses 
Mbh ~ 100 — 600 Mq (e.g., Heger & Woosley 2002). Such 
IMF would significantly reduce the mass ratio of merging 
BHs. In turn this would lower the value of gravitational wave 
recoil. 


4 DISCUSSION AND CONCLUSIONS 


and 

U|| = K cos (0 - 0o) ^ 5 (ai - qai) . (7) 

The fitting constants are A= 1.2x lO^kms”^, B= —0.93, H= 
(7.3±0.3) X 10®kms“^ and Kcos (0 — ©o) = (6.0±0.1) x 10^ , 
q is mass ratio of the merging BHs, Ui is the reduced spin pa¬ 
rameter and the orientation of the merger is specified by an¬ 
gles 0 and The authors show distribution of gravitational 
recoil for two models: when the spin parameters are chosen 
from a uniform distribution and when it is assumed that the 
BH spins are aligned with the orbital angular momentum. 
We accepted their model when BH spins are aligned with 
the orbital angular momentum (their fig. 2 , red region) and 
we used maximum of the kick velocity distribution (their fig. 
2 , red line) to assign values of kick velocities to the merging 
BHs in our model. 

Escape velocity is calculated as described in O’Leary & 
Loeb (2009). Circular velocity of a halo at the virial radius 
is (Barkana & Loeb 2001): 
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where Ac = IStt^ -|- 82d — 39d^, d = — 1 and = 

Hm(l + 2 )^/(Hm(l + 2 )^ + Ha) at the merger redshift. The 
authors assumed that the dark matter haloes are described 


In this paper we test under which conditions light BH seeds 
(100 Mq) placed in haloes of Millennium-H simulation 
(Boylan-Kolchin et al. 2009) and Millennium simulation 
(Springel et al. 2005) merger trees can grow to SMBHs with 
masses 10® Mq that have been observed at 2 ~ 7. 

We make merger trees which track dark matter halo 
merger history from redshift 2 = 23.79 to 2 = 6.2. We as¬ 
sume that each halo hosts one BH and that two BHs merge 
right after their host haloes merge. BH can grow in BH merg¬ 
ers and by gas accretion. We have distinguished between mi¬ 
nor and major mergers. In the case of a minor merger, mass 
of the newly formed BH is a simple sum of the previous 
BH masses. Major merger leads to the formation of a new 
BH and it triggers the gas accretion on to that BH. Final 
BH mass depends on the initial BH mass, Eddington ratio, 
radiative efficiency and quasar lifetime (equation 2.). Effec¬ 
tive radiative efficiency is fixed at e = 0.1 (Elvis, Risaliti & 
Zamorani 2002; Yu & Tremaine 2002; Davis & Laor 2011) 
and each accretion episode is limited to 50 Myr, which is 
~ Salpeter’s time for accretion at the Eddington limit. Ed¬ 
dington ratio is a free parameter in our model, but it has 
fixed value for each accretion episode in one simulation run. 
If a new major merger occurs before the maximum allowed 
time for accretion has passed, accretion is reset and the new 
accretion episode begins. Also, accretion can be stopped if 
BH mass exceeds 8 x 10“'* of the host halo mass, so we 
make sure that the gas reservoir never gets depleted. 

We combine Millennium-H and Millennium merger 
trees in order to have both: early halo formation history 
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with low-mass haloes (to track BH growth history); and a 
large box giving us abundant highest mass haloes later (in 
which 10® M© SMBH at z ~ 7 can be produced). First we 
place 100 Mq BH seeds in haloes of Millennium-II simula¬ 
tion which has 125 times better mass resolution to make 
BH growth history. Then we take most common BHs from 
Millennium-H simulation as seeds for Millennium simulation 
to produce 10® M© SMBH at 2 ~ 7. 

We run the set of semi-analytical simulations where we 
assign same initial masses of 100 Mq to all BH seeds and 
the same value for the effective Eddington ratio during the 
episodes of accretion. We investigate what value of the effec¬ 
tive Eddington ratio match two conditions: BH mass func¬ 
tion in our model cannot exceed BH mass function given by 
Willott et al. (2010b) and our merger trees need to produce 
10® Mq SMBH at 2 = 7. 

We find that remnants of Pop HI stars can produce 
10® M© SMBH at redshifts 2 = 7 if at each accretion episode 
they are able to accrete at the effective Eddington ratio of 
/fidd = 3.7. Recent observations have suggested that mod¬ 
erate super-Eddington accretion (1 < /Bdd < 10) might 

be possible (Kelly & Shen 2013; Du et al. 2014; Page et 
al. 2014; Novak 2013). 

The question is then how long can accretion maintain 
rates above the Eddington limit? Some authors argued that 
quasars lifetimes are much shorter then the Salpeter’s time- 
scale (e.g, Richstone et al. 1998; Wyithe & Padmanabhan 
2006). In that case, even massive BH seeds do not have 
enough time to grow to > 10® M© SMBH at redshift 2 > 7. 
A single BH seed with mass 10® M© requires nine e-folding 
times to produce > 10® M© SMBH at 2 = 7, while 100 M© 
seed requires 19 e-folding times for the same SMBH. Be¬ 
cause of this requirement some previous works rejected Pop 
HI star remnants as possible candidates for SMBH seeds 
(e.g., Johnson et al. 2012, 2013). 

In our model, growing SMBH at 2 = 7 does not require 
e-folding time larger than ~ 4 for 100 M© BH seeds and 
effective /Edd = 3.7. Thus, we have managed to reduce the 
time that BH needs to spend accreting. In our model, every 
major merger restarts the accretion, hence, instead of one 
BH constantly accreting for a long time (large number of 
e-folding times) we have several BHs in shorter (small num¬ 
ber of e-folding times) accretion episodes. This approach 
increases the importance of mergers on SMBH growth. 

Our model requires prolonged super-Eddington accre¬ 
tion which could produce strong feedback that stops the 
inflow of gas towards SMBH and disrupts the SMBH ac¬ 
cretion long before 50 Myr. This could occur if we would 
assume the classical ’thin disc’ accretion. However, we do 
not assume any accretion model in advance. The chosen val¬ 
ues for the radiative efficiency and the Eddington ratio in 
our model should be regarded as the effective values, aver¬ 
aged over 50 Myr. BH could also grow through a sequence 
of short-lived intermittent phases of super-Eddington accre¬ 
tion (/Edd > 3.7) without the feedback. Eddington ratio 
averaged over 50 Myr would then be /Edd = 3.7 and ef¬ 
fectively this model would produce the same SMBH as in 
the ’thin disc’ model. The absence of feedback in the ’slim 
disc’ model comes from fEdd = e ■ BH mass growth 
rate could stay the same as radiative efficiency decreases 
and /Edd increases. In principle, given a sufficiently low ef¬ 
ficiency, a super-Eddington BH may be emitting at sub- 


Eddington luminosity, thus eliminating BH feedback com¬ 
pletely (Volonteri & Silk 2015). 

In light of the recent observations of super- 
Eddington accretion mentioned above, we show 
that with a moderate super-Eddington accretion 
(/Edd = 3.7) averaged over 50 Myr, even low mass 
seeds (100 M©) could be progenitors of high redshift 
SMBHs (2 ~ 7). 

Note that our model is conservative for three reasons: 

(1) All BH seeds in our model have the same mass of 
100 M©. It is possible that Pop HI stars left BH seeds with 
masses up to 300 M© (Bromm et al. 2009), or even 1000 M© 
(Hirano et al. 2014). Larger mass of BH seeds would sig¬ 
nificantly reduce need for accretion which would, in turn, 
effectively lower the value of the effective Eddington ratio. 

(2) We combine two simulations which is also a con¬ 
servative approximation. We do not know the exact merger 
history of low mass haloes in Millennium simulation, but in¬ 
stead we take it from Millennium-H simulation. We always 
chose the most common BHs from Millennium-H simulation 
to be the seeds for newly formed haloes of Millennium Sim¬ 
ulation. If we would have chosen just a few BHs with mass 
greater than the most common mass, that would strongly 
affect BH growth history. This in turn would reduce the 
Eddington ratio necessary for SMBH growth. 

(3) In our model BH growth is limited by the amount 
of gas that BH can accrete (8 x 10“"^ of the host halo mass). 
This constraint comes from the BH mass scaling relations 
in the local Universe. However, the same scaling relation 
may not hold in the early Universe. Recently, Barnett et 
al. (2015) measured Mbh/Afbuige — 0.2 for quasar ULAS 
J1120-I-0641 at 2 = 7.1. They found that the BH was grow¬ 
ing much faster than the bulge relative to the mass ratio 
measured in the local Universe. If BHs at high redshift had 
more gas available for the accretion, BH growth would be 
much easier which would also reduce the value of the effec¬ 
tive Eddington ratio. 

Some authors have speculated that BH mergers have a 
limited role in SMBH growth due to gravitational wave recoil 
(e.g., Merritt et al. 2004; Volonteri 2007). We calculate BH 
kick velocities and compare it to the escape velocities of their 
host haloes. We find that BHs with low masses reside in low 
mass haloes and have their kick velocities larger than the 
escape velocities from their host haloes, so they are sensitive 
to the gravitational wave recoil. This is a consequence of a 
simple model where all BH seeds have the same mass which 
is highly unlikely. Using a proper initial BH mass function 
with a wide range of possible masses would lower the value 
of gravitational wave recoil. 

We note that, in order to get more accurate results, our 
model should be applied to the merger trees of a simulation 
that would have both resolution of Millennium-H simulation, 
and the box size as in Millennium simulation. 
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